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ABSTRACT 

A careful analysis of the HEAOl A2 2 — 10 keV full-sky map of the X-ray background (XRB) reveals 
clustering on the scale of several degrees. After removing the contribution due to beam smearing, the 
intrinsic clustering of the background is found to be consistent with an auto-correlation function of the 
form 3.6 ± 0.9 x 10~*9~^ where 9 is measured in degrees. If current AGN models of the hard XRB are 
reasonable and the cosmological constant-cold dark matter (ACDM) cosmology is correct, this clustering 
implies an X-ray bias factor of 6x ^ 2. Combined with the absence of a correlation between the XRB 
and the cosmic microwave background (CMB), this clustering can be used to limit the presence of an 
integrated Sachs- Wolfe (ISW) effect and thereby to constrain the value of the cosmological constant, 
f^A < 0.60 (95% C.L.). This constraint is inconsistent with much of the fl\ parameter space currently 
favored by other observations. Finally, we marginally detect the dipole moment of the diffuse XRB and 
find it to be consistent with the dipole due to our motion with respect to the mean rest frame of the 
XRB. The limit on the amplitude of any intrinsic dipole is Six /I < 5 x 10^"^ at the 95 % C.L. When 
compared to the local bulk velocity, this limit implies a constraint on the matter density of the universe 
of an° V&x(0) ;5 0.24. 

Subject headings: large-scale structure of the universe — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

The X-ray background (XRB) was discovered before the 
cosmic microwave background (CMB), but only now is its 
origin being fully understood. The hard (2 — 10 keV) 
XRB has been nearly completely resolved into individual 
sources; most of these are active galactic nuclei (AGN), 
but there is a minor contribution from the hot, intergalac- 
tic medium in rich clusters of galaxies (Rosati et al. 2002; 
Cowie et al. 2002; & Mushotzky et al. 2000). In ad- 
dition, the spectra of these faint X-ray sources are con- 
sistent with that of the "diffuse" XRB. If current models 
of the luminosity functions and evolution of these sources 
are reasonably correct, then the XRB arises from sources 
in the redshift range < z < 4, making them an impor- 
tant probe of density fluctuations intermediate between 
relatively nearby galaxy surveys {z <, 0.5) and the CMB 
(z ~ 1000). 

While there have been several attempts to measure large 
scale, correlated fluctuations in the hard XRB, these have 
only yielded upper limits or, at best, marginal detections 
(e.g. Barcons et al. 2000, Treyer et al. 1998 and refer- 
ences cited therein). On small scales, a recent correlation 
analysis of 159 sources in the Chandra Deep Field South 
survey detected significant correlations for separations out 
to 100 arcsec (Giacconi et al. 2001). (At the survey flux 
level, these sources comprise roughly two thirds of the hard 
XRB.) On much larger scales, a recent analysis by Scharf 
et al. (2000) claims a significant detection of large-scale 
harmonic structure in the XRB with spherical harmonic 
order 1 < ^ < 10 corresponding to structures on angular 
scales oi 9 ^ 10° The auto-correlation results we describe 



here complement this analysis, indicating clustering on an- 
gular scales of 3° to 10°, corresponding to harmonic order 
of ^ < 30. However, all three detections have relatively low 
signal to noise and require independent confirmation. 

The dipole moment of the XRB has received particular 
attention, primarily because of its relation to the dipole in 
the CMB, which is likely due to the Earth's motion with 
respect to the rest frame of the CMB. If this is the case, 
one expects a similar dipole in the XRB with an amplitude 
that is 3.4 times larger because of the difference in spectral 
indices of the two backgrounds (Boldt 1987). In the X-ray 
literature, this dipole is widely known as the Compton- 
Getting effect (Compton & Getting 1935). In addition, it 
is quite likely that the XRB has an intrinsic dipole due 
to the asymmetric distribution in the local matter density 
that is responsible for the Earth's peculiar motion in the 
first place. Searches for both these dipoles have concen- 
trated on the hard XRB, since at lower energies the X-ray 
sky is dominated by Galactic structure. There have been 
several tentative detections of the X-ray dipole (e.g. Scharf 
et al. (2000)), but these have large uncertainties. A firm 
detection of an intrinsic dipole or even an upper limit on 
its presence would provide an important constraint on the 
inhomogeneity of the local distribution of matter via a less 
often used tracer of mass and a concomitant constraint on 
cosmological models (e.g., Lahav, Piran & Treyer 1997). 

This paper is organized as follows. In section §2 we 
describe the hard X-ray map used in the analysis, the de- 
termination of its effective beam size, and cuts made to 
remove the foreground contaminants. In section §3, we 
describe the remaining large scale structures in the map 
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and the determination of their amphtudes. The dipole is 
of particular interest, and is the topic of section §4. The 
correlation function of the residual map and its implica- 
tions for intrinsic correlations are discussed in section §5. 
In section §6, we compare our results to previous obser- 
vations and discuss the cosmological implications of these 
results in §7. 

2. HEAOl A2 2- 10 KEV X-RAY MAP 

There has been much recent progress in understand- 
ing the X-ray background through instruments such as 
ROSAT, Chandra and XMM. However, these either have 
too low an energy threshold or have too small a field of 
view to study the large scale structure of the hard X- 
rays. The best observations relevant to large scale struc- 
ture are still those from the HEAOl A2 experiment that 
measured the surface brightness of the X-ray background 
in the 0.1 - 60 keV band (Boldt 1987). 

The HEAOl data set we consider was constructed from 
the output of two medium energy detectors (MED) with 
different fields of view (3° x 3° and 3° x 1.5°) and two high 
energy detectors (HED3) with these same fields of view. 
These data were collected during the six month period be- 
ginning on day 322 of 1977. Counts from the four detec- 
tors were combined and binned in 24,576 1.3° x 1.3° pixels. 
The pixelization we use is an equatorial quadrilateralized 
spherical cube projection on the sky, the same as used 
for the COBE satellite CMB maps (White and Stemwedel 
1992). The combined map has a spectral bandpass (quan- 
tum efficiency ^ 50%) of approximately 3—17 keV (Jahoda 
& Mushotzky 1989) and is shown in Galactic coordinates 
in Figure 0. For consistency with other work, all signals 
are converted to equivalent flux in the 2—10 keV band. 

Because of the ecliptic longitude scan pattern of the 
HEAO satellite, sky coverage and therefore photon shot 
noise are not uniform. However, the variance of the 
cleaned, corrected map, 2.1 x 10~^ {TOT counts s^^)^, 
is much larger than the variance of photon shot noise, 
0.8x 10-2 {TOT counts s'^f, where 1 TOT counts s'^ « 
2.1 X 10-"er5 s'^cm-"^ (Allen, Jahoda & Whitlock 1994). 
This implies that most of the variance in the X-ray map 
is due to "real" structure. For this reason and to reduce 
contamination from any systematics that might be corre- 
lated with the scan pattern, we chose to weight the pixels 
equally in this analysis. 

2.1. The point spread Junction 

To determine the level of intrinsic correlations, we must 
account for the effects of beam smearing and so it is es- 
sential to characterize the point spread function (PSF) of 
the above map. The PSF varies somewhat with position 
on the sky because of the pixelization and the asymmet- 
ric beam combined with the HEAOl scan pattern. We 
obtained a mean PSF by averaging the individual PSFs 
of 60 strong HEAOl point sources (Piccinotti 1982) that 
were located more than 20° from the Galactic plane. The 
latter condition was imposed to avoid crowding and to ap- 
proximate the windowing of the subsequent analysis (see 
§2.2). The composite PSF, shown in Figure ||, is well fit 
by a Gaussian with a full width, half maximum (FWHM) 
of 3.04°. 

As a check of this PSF, we generated Monte Carlo maps 
of sources observed with 3° x 3° and 3° x 1.5° (FWHM) 



triangular beams appropriate for the A2 detectors (Shafer 
1983) and then combined the maps with quadcubed pix- 
elization as above. The resulting average PSF from these 
trials is also well fit by a Gaussian with a FWHM of 2.91°, 
i.e., about 4.5% less than that in the Figure Consider- 
ing that the widths of the triangular beams given above 
are nominal, that the triangular beam pattern is only ap- 
proximate (especially at higher energies), and that we did 
not take into account the slight smearing in the satellite 
scan direction (Shafer 1983), the agreement is remarkably 
good. In the following analysis, we use the 3.04° fit de- 
rived from the observed map; however, changing the PSF 
FWHM by a few percent does not significantly affect the 
results of this paper. 

2.2. Cleaning the map 

To remove the effects of the Galaxy and strong extra- 
galactic point sources, some regions of the map were ex- 
cluded from the analysis. The dominant feature in the 
HEAO map is the Galaxy (see Figure P so all data within 
20° of the Galactic plane or within 30° of the Galactic 
center were cut from the map. In addition, large regions 
(6.5° x 6.5°) centered on 92 discrete X-ray sources with 
2—10 keV fluxes larger than 3 x IQ-^^erg s'^cm"^ (Pic- 
cinotti 1982) were removed from the maps. Around the 
sixteen brightest of these sources (with fluxes larger than 

1 X lO^^^ergi s~^cm.~^) the cut regions were enlarged to 
9° X 9°. Further enlarging the area of the excised regions 
had a negligible effect on the following analysis so we con- 
clude that the sources have been effectively removed. The 
resulting "cleaned" map (designated Map A) has a sky 
coverage of 55.5% and is our baseline map for further cuts. 

To test the possibility of further point source contami- 
nation, we also used the ROSAT All-Sky Survey (RASS) 
Bright Source Catalog (Voges et al. 1996) to identify rel- 
atively bright sources. While the RASS survey has some- 
what less than full sky coverage (92%), it has a relatively 
low flux limit that corresponds to a 2 — 10 keV flux of 
~ 2 X lO"^"^ erg s^^ cm^^ for a photon spectral index of 
a = —2. Every source in the RASS catalog was assigned a 

2 — 10 keV flux from its B-band flux by assuming a spectral 
index of— 3<q;<— las deduced from its HR2 hardness 
ratio. For fainter sources, the computed value of a is quite 
uncertain; if it fell outside the typical range of most X-ray 
sources, — 3 < a < — 1, then a was simply forced to be 
— 1 or —3. It is clear that extrapolating RASS flux to the 
2 — 10 keV band is not accurate, so one must consider the 
level to which sources are masked with due caution. How- 
ever, we are only using these fluxes to mask bright sources 
and so this procedure is unlikely to bias the results. 

We considered maps where the ROSAT sources were 
removed at three different inferred 2—10 keV fiux thresh- 
olds. First, we identified sources with fluxes exceeding the 
Piccinotti level, 3 x l{)~^^erg s~^cm~^ . Thirty-four addi- 
tional, high Galactic latitude RASS sources were removed, 
resulting in a map with sky coverage of 52% (designated 
Map B). In order to compare more directly with the re- 
sults of Scharf et al. (2000) (see §6) we removed sources at 
their flux level, 2 x \Q~^^erg s~^cm~^. The map masked 
in this way has 47% sky coverage (compared to the 48% 
coverage of the Scharf et al. analysis) and is designated 
Map C. Finally, to check how sensitive our dipole results 
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are to the particular masking of the map, we lowered the 
flux cut level to 1 x IQ'^^erg s~^cm^^, which reduced the 
sky coverage to 34%. The map resulting from this cut is 
designated Map D in Table 3. 

As an alternative to using the RASS sources, the map 
itself was searched for "sources" that exceeded the nearby 
background by a specified amount. Since the quad-cubed 
format lays out the pixels on an approximately square 
array, we averaged each pixel with its eight neighbors 
and then compared this value with the median value of 
the next nearest sixteen pixels (ignoring pixels within the 
masked regions). If the average flux associated with a 
given pixel exceeded the median flux of the background 
by a prescribed threshold, then all 25 pixels (6.5° x 6.5°) 
were removed from further consideration. For a thresh- 
old corresponding to 2.2 times the mean shot noise in the 
map approximately 120 more "sources" were identified and 
masked resulting in a sky coverage of 42%. This map is 
labeled Map E in Table 3. Finally, we used an even more 
aggressive cut corresponding to 1.75 times the mean shot 
noise which resulted in a masked map with 33% sky cov- 
erage. This map is labeled Map F. 

3. MODELING THE LOCAL LARGE-SCALE STRUCTURE 
3.1. Sources of large scale structure 

There are several local sources of large-scale structure 
in the HEAO map which can not be eliminated by mask- 
ing isolated regions. These include diffuse emission from 
the Galaxy, emission (diffuse and/or faint point sources) 
from the Local Supercluster, the Compton-Getting dipole, 
and a linear time drift in detector sensitivity. Since none 
of these are known a priori, we fit an eight parameter 
model to the data. Of course, the Compton-Getting dipole 
is known in principle if one assumes the kinetic origin of 
the dipole in the cosmic microwave background; however, 
there may also be an intrinsic X-ray dipole that is not 
accounted for. (See §4 below.) Only one correction was 
made a priori to the map and that was for the dipole 
due to the Earth's motion around the sun; however, this 
correction has a negligible effect on the results. A more 
detailed account of the model is given in Boughn (1999). 

The X-ray background has a diffuse (or unresolved) 
Galactic component which varies strongly with Galactic 
latitude (Iwan et al. 1982). This emission is still signif- 
icant at high Galactic latitude {bu > 20°) and extrap- 
olates to ~ 1% at the Galactic poles. We modeled this 
emission in two ways. The first model consisted of a lin- 
ear combination of a secant law Galaxy with the Haslam 
408 GHz full sky map (Haslam et al. 1982). The lat- 
ter was included to take into account X-rays generated by 
inverse Compton scattering of CMB photons from high 
energy electrons in the Galactic halo, the source of much 
of the synchrotron emission in the Haslam map. As an 
alternative Galaxy model we also considered the two disk, 
exponentially truncated model of Iwan et al. (1982). Our 
results are independent of which model is used. 

In addition to the Galactic component, evidence has 
been found for faint X-ray emission from plane of the Lo- 
cal Supercluster (Jahoda 1993, Boughn 1999). Because 
of its faintness, very detailed models of this emission are 
not particularly useful. The model we use here is a simple 
"pillbox" , i.e uniform X-ray emissivity within a circular 
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disk of thickness equal to a 1/4 of the radius and with 
its center located 4/5 of a radius from us in the direction 
of the Virgo cluster (see Boughn 1999 for details). The 
amplitude of this emission, while significant, is largely in- 
dependent of the details of the model and, in any case, has 
only a small effect on the results. 

Time drifts in the detector sensitivity can also lead to 
apparent structure in the reconstructed X-ray map. At 
least one of the A-2 detectors changed sensitivity by ~ 1% 
in the six month interval of the current data set (Jahoda 
1993). Because of the ecliptic scan pattern of the HEAO 
satellite, this results in a large-scale pattern in the sky 
which varies with ecliptic longitude with a period of 180°. 
If the drift is assumed to be linear, the form of the resulting 
large-scale structure in the map is completely determined. 
A linear drift of unknown amplitude is taken into account 
by constructing a sky map with the appropriate structure 
and then fitting for the amplitude simultaneously with the 
other parameters. We investigated the possibility of non- 
linear drift by considering quadratic and cubic terms as 
well; however, this did not significantly reduce the of 
the fit nor change the subsequent results. 

3.2. Modeling the maps 

The eight parameters that characterize the amplitude 
of these structures are used to model the large-scale struc- 
ture in the HEAO map. Let the X-ray intensity map be 
denoted by the vector I, where the element It is the inten- 
sity in the i'^ pixel. The observed intensity is modelled as 
the sum of eight templates with amplitudes described by 
the eight dimensional vector a, 

I = Xa + n (3-1) 

where X is an 8 x ripi^ matrix whose elements are the val- 
ues of each template function at each pixel of the map. 
As discussed above, these template functions include: a 
uniform map to represent the monopole of the X-ray back- 
ground; the three components of a dipole (in equatorial co- 
ordinates); the large-scale pattern resulting from a linear 
instrumental gain drift; a Galactic secant law; the Haslam 
408 GHz map ; and the amplitude of the "pillbox" model 
of the local supercluster. The noise vector n is assumed 
to be Gaussian distributed with correlations described by 

As discussed above (§2), we chose to weight each pixel 
equally since the shot noise is considerably less than the 
"real" fluctuations in the sky. For the purposes of fitting 
the map to the above model we consider both photon shot 
noise and fluctuations in the XRB (see Figure H) to be 
"noise" . This noise is correlated and a minimum fit 
must take such correlations into account. However, for 
simplicity, we ignore these correlations when finding the 
best fit model amplitudes and perform a standard least 
squares fit by minimizing 1 1 - Xap on the cleaned HEAO 
map. From the standard equations of linear regression the 
values of the parameters that minimize this sum are 

a = B-^X^I (3-2) 

where B — X^X is a symmetric eight by eight matrix. 
This would be the maximum likelihood estimator if the 
correlation matrix were uniform and diagonal. Though 
this fit ignores correlations in the errors, it is unbiased 
and is likely to be very close to the minimum (maxi- 
mum likelihood) fit, since the noise correlations are on a 
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much smaller scale than the features we are attempting to 
fit. 

The correlated nature of the noise cannot be ignored 
when computing the uncertainties in the fit since there 
are far fewer noise independent data points than there are 
pixels in the map. It is straightforward to show that errors 
in the estimated parameters 5a are given by 

(SaSa^) = B-^X'^CXB-^ (3-3) 

This error is likely to be only slightly larger than would 
be the case for the maximum likelihood estimator. C is a 
combination of the uncorrelated shot noise and the corre- 
lated fluctuations indicated in Figure ^. We assume it to 
be homogeneous and isotropic, i.e., that Cij depends only 
on the angular separation of the i and j pixels. 

Table lists the values and errors of the parameters 
fit to Map C (see §2). Instrument time drift, the Galaxy 
and structure associated with the local supercluster all ap- 
pear to be significant detections. The dipole is detected at 
about the 2 a level and is consistent with that expected 
for the Compton-Getting dipole (see Table Table |^ 
lists the elements of the normalized correlation matrix of 
the fit parameters and it is apparent that the parameters 
are largely uncorrelated. This was supported by fits that 
excluded some of the parameters (see §4) . 

To compute the true = (I-Xa)^C'-i(I-Xa) of the 
fit requires inverting C which is a 11, 531 x 11, 531 matrix. 
Instead we compute an effective reduced using 

xlff^^iI-XafD-\l-Xa) (3-4) 

where D is the diagonal part of the correlation matrix. 
Da = cr"^ i + cr^, (Ts,i is the shot noise in the i*'' pixel, crf,^ 
is the variance of the fluctuations in the XRB, and N is 
the number of pixels minus eight, the number of degrees 
of freedom in the fit. The shot noise in a given pixel is 
inversely proportional to the number of photons received 
and we assume is inversely proportional to the coverage of 
that pixel. This is approximately true since all the non- 
fiagged pixels are exposed to approximately the same fiux. 
We find Xe//^ = 1-00 for this fit, which we take as an indi- 
cation that we have properly characterized the amplitude 
of the noise so that the errors quoted in the table have nei- 
ther been underestimated nor overestimated. However, it 
should be emphasized that Xe//^ is not to be interpreted 
statistically as being derived from a x^ distribution. 

The residual maps show very little evidence for struc- 
ture on angular scales 9 > 10° above the level of the noise, 
{5P)/P ~ 10~^ wher_e SI are the residual fiuctuations in 
X-ray intensity and / is the mean intensity (see Figure 
H). Since all the components of the model have significant 
structure on large angular scales, it appears that these 
particular systematics have been effectively eliminated. 

4. THE DIPOLE OF THE X-RAY BACKGROUND 

The dipole fit to the map is consistent with the 
Compton-Getting dipole and there is no evidence for any 
additional intrinsic dipole in the XRB. To make this more 
quantitative, we corrected the maps for the predicted 
Compton-Getting dipole and fit the corrected map for any 
residual, intrinsic dipole. These dipole fit parameters are 
also included in Table |^. 



Leaving out any individual model component, such as 
the time drift, the galaxy, Haslam or supercluster tem- 
plate, made little difference in the amplitude of the fit 
dipole. This is, perhaps, not too surprising since the 
Galaxy and time drift models are primarily quadrupolar 
in nature and the pancake model, while possessing a sig- 
nificant dipole moment, has a relatively small amplitude. 
All such fits were consistent with the Compton-Getting 
dipole alone. Even when all four of these parameters were 
excluded from the fit, the dipole amplitude increased by 
only 0.004 TOT counts with a direction that was was 
33° from that of the CMB dipole. The effective x^ for the 
four parameter fit was, however, significantly worse, i.e., 

Xe//2 = 1.05. 

Table g lists the amplitude and direction of the dipole 
fit to Map C along with the fits to Maps D, E, and F. 
All of these fits are consistent with amplitude and direc- 
tion of the Compton-Getting dipole (as inferred from the 
CMB dipole) which is also indicated in the Table. The 
effective x^s of these fits range from 0.99 to 1.01, again 
indicating that the amplitude of the noise is reasonably 
well characterized. No errors are given for these quanti- 
ties for reasons that will be discussed below. In order to 
check for unknown systematics, we performed dipole fits 
to a variety of other masked maps with larger Galaxy cuts 
as well as cuts of the brighter galaxies in the TuUy Nearby 
Bright Galaxy Atlas (TuUy 1988). The details of these 
cuts are discussed in Boughn (1999); however, none had a 
significantly different dipole fit. 

Since all of these dipoles are consistent with the 
Compton-Getting dipole we also fit these maps with a six 
parameter fit in which the dipole direction was constrained 
to be the direction of the CMB dipole. The dipole ampli- 
tude of these fits and errors computed according to Eq. 
(3-3) are also given in Table |3[ 

Even though we find no evidence for an intrinsic dipole 
in the XRB, it would be useful to place an upper limit 
on its amplitude. We define the dimensionless dipole by 
writing the first two moments of the X-ray intensity as 

/(fi) = /(l + A-n), (4-1) 

where A is a vector in the direction of the dipole. There 
are various approaches one could take to find an upper 
limit, and the problem is complicated somewhat because 
the error bars are anisotropic (see Table 1). The dipole 
in the y direction is less constrained than in the other di- 
rections because of the anisotropic masking of the map. 
Here we take the limits on the individual components of 
the intrinsic dipole and marginalize over the dipole direc- 
tion to obtain a distribution for its amplitude. For this, 
we use a Bayesian formalism and assume a uniform prior 
on the amplitude, |A|. We find A < 0.0052 at the 95 % 
C.L. If the direction of the dipole is fixed to be that of the 
CMB dipole then the 95% C.L. upper limits on the dipole 
amplitudes fall in the range 0.0030 to 0.0043 for the fits 
listed in Table |. 

The same sort of problem arises when trying to attach 
an error bar to the amplitude of the dipole fits to the maps 
which include the Compton-Getting dipole. It seems clear 
from Table |l] that we find evidence for a dipole at the 2 
(7 level. This is supported by the six parameter fits of Ta- 
ble 1^, where the various maps indicate positive detections 
at a 2 to 3 cr level. However, in the eight parameter fits, 
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Parameter 


Uncorrected 


Corrected for C-G 


ai 


background 


oZa.b ± i.y 


same 


(12 


X dipole 


-i.l I ± U.Dz 


-U.z4 ± U.DZ 




y dipole 


-0.38 ± 0.98 


-0.68 ± 0.98 


ai 


z dipole 


-0.52 ± 0.69 


-0.34 ± 0.69 


05 


time drift 


7.15 ± 1.23 


same 


oe 


secant law 


3.28 ± 0.84 


same 


a^ 


Haslam map 


0.03 ± 0.08 


same 


as 


Supercluster 


4.11 ± 1.35 


same 



Table 1 



Eight fit parameters for Map C 
0.01 TOT count s-^(4.5 deg'^)~^ ~ 1.54 x 10~'''erg s 
map corrected for the Compton-Getting (C-G) dipole. 



(sources brighter than 2 x 10 erg s cm removed). The units are 
^cm~^. Fits are shown both for the original map and for the 





ai 


0-2 


03 


04 


ag 


ae 


ar 


as 


ai 


1.0 


0.0 


-0.5 


0.1 


-0.4 


-0.6 


-0.1 


-0.3 




0.0 


1.0 


-0.1 


0.2 


0.3 


-0.1 


0.0 


0.1 


03 


-0.5 


-0.1 


1.0 


0.1 


0.0 


0.7 


0.0 


-0.3 


04 


0.1 


0.2 


0.1 


1.0 


-0.1 


0.0 


0.0 


-0.1 


05 


-0.4 


0.3 


0.0 


-0.1 


1.0 


0.0 


0.0 


0.2 


ae 


-0.6 


-0.1 


0.7 


0.0 


0.0 


1.0 


0.0 


-0.4 


ar 


-0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


-0.3 


0-8 


-0.3 


0.1 


-0.3 


-0.1 


0.2 


-0.4 


-0.3 


1.0 



Table 2 



Normalized correlation coefficients for the fit parameters in Table |l|. 



Map 


Dipole 


hi 


bii 


Constrained Dipole 


Map C 


0.0133 


309" 


39° 


0.0117 ± 0.0064 


Map D 


0.0218 


300° 


33° 


0.0184 ± 0.0062 


Map E 


0.0150 


296° 


50° 


0.0148 ± 0.0059 


Map F 


0.0190 


283° 


44° 


0.0184 ± 0.0064 


C-G 


0.0145 


264° 


48° 


0.0145 



Table 3 



The dipole amplitude and directions are from the 8-parameter fits in TOT units and Galactic coordinates. Map C is for 
the map of Table |^; Map D is masked at a source level of ^ 1 x lO^^^erg s^^cm^^; Map E is masked with internal source 
identification; and Map F is masked with a lower level of internal source identification (see §2 for full details). Also listed 
are the amplitude and direction of the Compton-Getting dipole (G-P) as inferred from the CMB dipole. The constrained 
amplitudes are for dipole models fixed to the direction of the CMB dipole. 



the dipole amplitude is a non-linear combination of the 
fit components. There are two approaches we can take in 
converting the three dimensional limits to a limit on the 
dipole amplitude. We can either fix the direction in the di- 
rection of the CMB dipole, which results in the constrained 
limits shown in Table pi Alternatively, we can marginalize 
over the possible directions of the dipole, which will nec- 
essarily result in weaker limits than when the direction is 
fixed. This is particularly true here, where the direction of 
greatest uncertainty in the dipole measurement is roughly 
orthogonal to the expected dipole direction. In the case of 
the Compton-Getting dipole, there is a strong prior that 
it should be in the CMB dipole direction, so our limit is 
stronger than it would be if we did not have information 
about the CMB. 



In addition to the upper limit on the intrinsic dipole 
amplitude, we can also constrain the underlying dipole 
variance, which can, in turn, be used to test theoretically 
predicted power spectra. While the observed amplitude is 
related to the dipole variance, (A^) = 3a\, there is large 
uncertainty due to cosmic variance. The dipole represents 
only three independent samplings of a a- To constrain a a, 
we again take a Bayesian approach and calculate the like- 
lihood of observing the data given the noise and a a, 

7'(A|aA) oc n e-'^^/^^-.'+^A) (^2 ^ ^|)-i/2^ (4.3) 

where the product is over the three spatial directions 
and we have ignored the small off-diagonal noise corre- 
lations (see Table ^). With a uniform prior on a a, its 
posterior distribution implies a 95% C.L. upper limit of 
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(TA < 0.0064. This is twice as high as would be inferred 
from the Umit on the dipole because of the significant tail 
in the distribution due to cosmic variance. The limit im- 
phed by the dipole (of Map C), cta = A/VS < 0.0030 is 
at the 80% C.L. The difference between the limits arises 
because occasionally a small dipole can occur even when 
the variance is large. 

The bottom line is that we have detected the dipole in 
the XRB at about the 2 a level and that it is consistent 
with the Compton-Getting dipole. There is no evidence 
for any other intrinsic dipole at this same level. We will 
discuss the apparent detection of an intrinsic dipole by 
Scharf et al. (2000) in §6. 

5. CORRELATIONS IN THE X-RAY BACKGROUND 

A standard way to detect the clustering of sources (or 
of the emission of these sources) is to compute the auto- 
correlation function (ACF), defined by 

^{e)^^Y.^h-I){lj-I)INe (5-1) 

where the sum is over all pairs of pixels, i, j, separated by 
an angle 6, Ii is the intensity of the ith pixel, / is the mean 
intensity, and Ng is the number of pairs of pixels separated 
by 0. Figure || shows the ACF of the residual map after 
being corrected with the 8-parameter fit and for photon 
shot noise in the 9 = 0° bin. The error bars are highly 
correlated and were determined from Monte Carlo trials 
in which the pixel intensity distribution was assumed to 
be Gaussian with the same ACF as in the figure. There is 
essentially no significant structure for 9 > 13° once local 
structures have been removed, as is evident in Figures |[ 
I and I 

It is clear from Figure Q that the residuals of Map A 
possess significant correlated structure. It must be deter- 
mined how much, if any, is due to clustering in the XRB 
and how much is simply due to smearing by the PSF of the 
map. It is straightforward to show that an uncorrelated 
signal smeared by a Gaussian PSF, PSF{9) oc /2<^p^ 

results in an ACF of the form u;{0) cx e"'^'/"'^' where 
(Tp = 1.29° is the Gaussian width of the PSF in Figure 
H i^FWHM — ScTp In 2.) The dashed curve in Figure^ is 
essentially this functional form, modified slightly to take 
into account the pixelization. In the plot, its amplitude 
has been forced to agree with 6* = 0° data point, while a 
maximum likelihood fit results in an amplitude about 5% 
lower (a consequence of the correlated noise). For 9 > 3° , 
the ACF of the data clearly exceeds that accountable by 
beam smearing and this excess is even more pronounced 
with the maximum likelihood fit. The reduced for the 
fit to the first eight data points {9 ^ 9°) is = 18.6 for 
six degrees of freedom which is another measure of the ex- 
cess structure between 3 and 9 degrees. Note, this is a two 
parameter fit since the photon shot noise (which occurs 
only a,t 9 = 0°) is also one of the parameters. 

While it is apparent that there is some intrinsic corre- 
lation in the X-ray background, it cannot be estimated 
by the residual to the above two parameter fit since that 
overestimates the contribution of beam smearing in or- 
der to minimize x^. Instead, we also include in the fit a 
form for the intrinsic correlation and find its amplitude as 
well. Since the signal to noise is too small to allow a de- 
tailed model of the intrinsic clustering, we chose to model 



it with a simple power-law uj{9) = (9q/9Y. This form 
provides an acceptable fit to the ACF of both radio and 
X-ray surveys on somewhat smaller angular scales (Cress 
& Kamionkowski 1998, Soltan et al. 1996, Giacconi et al. 
2001). This intrinsic correlation was then convolved with 
the PSF and applied to the quadcube pixelization of the 
map (e.g., Boughn 1998). 

Finally, it is important to take into account the effects of 
the 8-parameter fit used to remove the large-scale structure 
as discussed in §3. If the X-ray background has intrinsic 
structure on the scale of many degrees, the 8-parameter fit 
will tend to remove it in order to minimize x^- Since the 
model is composed of relatively large scale features, the 
greatest effect is expected for the largest angles. The sig- 
nificance of this effect was determined by generating Monte 
Carlo trials assuming a Gaussian pixel intensity distribu- 
tion with the same ACF as in Figure ||. The 8-parameter 
model was then fit and each trial map was corrected ac- 
cordingly. The ACFs computed for these corrected maps 
indicate that, as expected, the value of the ACF is signifi- 
cantly attenuated for larger angles. The attenuation factor 
for 6* = 9° is already 0.55 and decreases rapidly for larger 
angles. The errors indicated in Figure |^ were also deter- 
mined from these Monte Carlo trials and, as mentioned 
above, are highly correlated. 

We model the auto-correlation function as a sum of three 
templates and fit for their best amplitude. Much of the 
analysis parallels the discussion for the fits of large scale 
structure in the map, with the exception that here the 
number of bins is small enough that it is simple to cal- 
culate the maximum likelihood fit. Again, we model the 
observed correlation function vector as = Wc + n^, 
where is a 3 x n^in matrix containing the templates 
for shot noise (ws), beam smearing {lupsf), and intrinsic 
correlations in the XRB {comtr)- The amplitudes are given 
by the three element vector c and the noise is described 
by the correlation matrix determined from Monte Carlo 
trials, = (n^nj). 

Shot noise contributes only to the first bin (zero sepa- 
ration) of the ACF and has amplitude given by its vari- 
ance, ujs- Beam smearing contributes to the ACF with a 
template that looks like the beam convolved with itself, 
so appears as a Gaussian with a FWHM a factor of \/2 
larger than that of the beam and has amplitude denoted 
as LOpsf. Finally, the intrinsic correlations are modelled 
as {9o/9Y, which is then smoothed appropriately by the 
beam. Its amplitude is denoted by its inferred correlation 
at zero separation, LOintr- We fit using a range of indices, 
0.8 > e > 1.6, which cover the range of theoretical models 
of the intrinsic correlation. Both the PSF template and 
the intrinsic template are modified to include the effects 
of the attenuation at large angles as discussed above. 

Minimizing with respect to the three fit parameters, 
c, results in the maximum likelihood fit to the model if 
one assumes Gaussian statistics. This assumption is rea- 
sonable by virtue of the central limit theorem since each 
data point consists of the combination of the signals from 
a great many pixels, each of which is approximately Gaus- 
sian distributed. In the presence of correlated noise, is 
defined by 

X^^{u;-WcfC-\Lo~Wc) (5-2) 
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It is straightforward to show that the value of the param- 
eters that minimize £^re given by 

c = n-^W^C-^oj (5-3) 

where = C~^W . Because of the large attenuation of 
the ACF at large angles, we chose to fit to only those data 
points with 0i < 9°, i.e., i < 8, even though there appears 
to be statistically significant structure out to 6* ^ 13°. The 
results of the fit of the e = 1 model to the ACF of the resid- 
uals of Map A are listed in Table ^ and plotted in Figure 

I 

It is also straightforward to show that the correlation 
matrix of the fit parameters is given by, {6cnScm) = 
so the errors given in Table ^ are given by cr^^^ = 
and the normalized correlation coefficients by r„m = 
f2~4/(r2~^ri~^)^/^. The correlations are as expected, i.e., 
LOintr and ujpsf are highly correlated, while lu^ is relatively 
uncorrelated with the other two parameters. 

From the results in Table ^ it appears that intrin- 
sic correlations in the X-ray background are detected 
at the 4 a level for sources with flux levels below 3 x 
10~^^erg s~^cm~'^. Of course, if we have not success- 
fully eliminated sources with fluxes larger than this, then 
the clustering amplitude might well be artificially inflated. 
It was for this reason that we used the RASS catalog 
to identify and remove additional sources with intensities 
^ 3 X lO^^^erg s^^cm^'^, the result of which was Map B 
(see §2). The fits to Map B are also listed in Table |. The 
clustering amplitude, comtr, of the fits to this modified map 
were only 11% less than that of Map A, i.e., considerably 
less than 1 a. The x^s of both fits are acceptable. 

These results are not very sensitive to the attenuation 
corrections. If they are removed from the model, the re- 
sulting amplitude of the fit clustering coefficient only de- 
creases by ~ 20%, less than 1 a. It should be noted that 
the corrections were relatively small (an average attenua- 
tion factor of 0.83 ranging from 1.0 at 6* = 0° to 0.55 for 
9 = 9°) and in all cases, the attenuation was less than the 
error bar of the corresponding data point. If data points 
with 9 > 9° are included, the fits become more sensitive 
to the attenuation corrections which are, in turn, quite 
sensitive to the 8-parameter fit. 

Figure ^ is a plot of the residuals of the fit to the ACF 
of Map A for 6* < 9° and of the uncorrected ACF from 
10° to 180°. The vertical scale is the same as for Figure 
1^. The rms of these 140 data points is 1 x 10~^ and it is 
clear that there is very little residual structure at levels ex- 
ceeding this value. The observed correlation function for 
9 > 10° in Figure || is entirely consistent with the noise 
levels determined from the Monte Carlo simulations: the 
rms oiuj/a is 1.03, indicating that there is no evidence for 
intrinsic fluctuations on these scales. We also take this as 
an indication that the errors are reasonably well character- 
ized by the Monte Carlo calculation. The variance of the 
photon shot noise, CTs^ = oJsP , is consistent with that ex- 
pected from photon counting statistics only (Jahoda 2001). 
It should be noted that not only the shape but the am- 
plitude of the beam smearing contribution, ujpsf, can be 
computed from source counts as a function of 2 — 10 keV 
flux. We will argue in §6 that, while our fltted value is 
consistent with current number counts, the latter are not 
yet accurate enough to correct the data. 

Figure || shows the model of the intrinsic clustering, 



^intr, compared to the data with both the shot noise and 
PSF component removed. The model curve is not plotted 
beyond 9 = 9° since the attenuation factor due to the 8- 
parameter flt corrections are large and uncertain at larger 
angles. The amplitude of ujintr is sensitive to the exponent 
in the assumed power law for the intrinsic correlations. For 
example, the fit amplitude for a 9^^'^ power law the am- 
plitude is a factor of ~ 2 larger than for a 6~^-^ power law. 
However, for a range of fits with 0.8 < e < 1.6, the values 
of ujintr{9) a,t 6 = 4.5° are all within ±3% of each other. 
Therefore, we chose to normalize the X-ray ACF at 4.5° 
when comparing to cosmological models (see §7.1). The 
X^s are reasonable for all these fits. 

The bottom line is that there is fairly strong evidence 
for intrinsic clustering on these angular scales at the level 
oi LUintr ^ 3.6 X 10~'^9~^ (see §6.1). While the exponent 
is e ~ 1, it is not strongly constrained. The implications 
of the intrinsic clustering of the X-ray background will be 
discussed in §7. 

6. COMPARISONS WITH PREVIOUS WORK 

6.1. Clustering 

As mentioned above, the component of the ACF due to 
beam smearing, lopsf, can be determined with no free pa- 
rameters if the fiux limited number counts of X-ray sources 
are known. While such counts are still relatively inaccu- 
rate for our purposes, we did check to see if our results are 
consistent with current data. Over a restricted flux range, 
the X-ray number counts, A^(< S), are reasonably approx- 
imated by a power law, i.e., A^(< S) = K where S is 
the flux of the source. It is straightforward to show that 
the variance of flux due to a Poisson distribution of sources 
is 

45f(0) = ^^p' A^^K S^-y{2 - 7) (6-1) 

where S is the upper limit of source flux, ap is the Gaussian 
width of the PSF and A^^ is the flux of a point source that 
results in a peak signal of 1 TOT count in our compos- 
ite map. In our case, A = 2.05 x lO^^erg^^ s cm?. Using 
the BeppoSAX 2 — 10 kev number count data of Giommi, 
Perri, & Fiore (2000), the Chandra data of Mushotzky 
et al. (2000), and the HEAOl A2 data of Piccinotti et 
al. (1982), we constructed a piecewise power-law N{< S) 
for the range 6 x 10^^^ < S" < 3 x 10~^^er5' s~^cm^^ and 
computed iOpspiO) = (^^psp/P to be 8.4 x 10"'*. The close 
agreement of this value with those in Table ^ is fortuitous 
given that the value depends most sensitively on the num- 
ber counts at large fluxes which are the most unreliable, 
typically accurate only to within a factor of two. However, 
it is clear that the lopsf of Table ^ are quite consistent 
with the existing number count data. 

The results of §5 clearly indicate the presence of intrin- 
sic clustering in the XRB. Mindful that the detection is 
only ~ 4 CT, we tentatively assume that the ACF has an 
amplitude of ujintr ^ 2.5 x 10~^ (the average of the two 
values in Table |^) and is consistent with a, 9~^ functional 
dependence. It is straight forward to relate this to the 
underlying correlation amplitude, 9f). The variance of cor- 
relations smoothed by a beam of Gaussian size a is given 

by 

r(l-e/2) (9^y 
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Map 


Us X 10'' 


LOpsF ^ 


Wmtr X lO'* 


xi 


ri2 


ri3 


r23 


A 


7.71 ± 0.15 


8.63 ± 0.59 


2.64 ± 0.65 


2.3/5 


-0.19 


-0.04 


-0.83 


B 


7.57 ± 0.15 


8.22 ± 0.59 


2.33 ± 0.65 


3.6/5 


-0.19 


-0.04 


-0.83 



Table 4 



Fit model parameters for Map A with 55% sky coverage and Map B with 52% sky coverage after removing additional 
ROSAT sources (see §2). The intrinsic fluctuations are modelled a,s uj (x . 



When e = 1, then Oq ~ 2(7UJintr /t^^^"^ ■ Using this, we find 

cjxflsW - 3.6 X 10"'' 61-1 (6-3) 

where 6 is measured in degrees and the normalization is 
such that wxflB(O) = {Sl'^)/P. 

For comparison, this amplitude is about a factor of three 
below the 2 a upper limit determined by Carrera et al. 
(1993) obtained with Ginga data (4 — 12 keV) for angu- 
lar scales between 0.2 and 2.0 degrees. The detection of 
a significant correlation in the HEAOl A2 data at the 
level of 3 X 10"^ at 6* = 10° by Mushotzky and Jahoda 
(1992) was later attributed to structure near the super- 
Galactic plane (Jahoda 1993). To check for this effect in 
the present analysis, we masked all pixels within 15 and 20 
degrees from the super-Galactic plane. The results were 
indistinguishable from those of Table 0. It is interesting 
that the clustering indicated in Eq. (6-3) is consistent with 
this level of fluctuations at 10°; however, our sensitivity 
has begun to decline significantly at 10° due to the fit for 
large scale structures. A correlation analysis of ROSAT 
soft X-ray background by Soltan et al. (1996) detected 
correlations about an order of magnitude larger than in- 
dicated in Eq. (6-3). Even considering that the ROSAT 
band (0.5 - 2.0 keV) is distinct from the HEAO band, it 
is difficult to imagine that the two correlation functions 
could be so disparate unless the lower energy analysis is 
contaminated by the Galaxy. 

While there has yet to be a definitive detection of the 
clustering of hard X-ray sources, a recent deep Chandra 
survey of 159 sources shows a positive correlation of source 
number counts on angular scales of 5 to 100 arcsec (Giac- 
coni et al. 2001). Although the signal to noise is low and 
dependent on source flux, the implied number count ACF 
is roughly consistent with ujN{d) '--^ 3 x 10~^d~^ where 
ujn{0) = {5N'^)/N'^ and N is the mean surface density 
of sources. This is consistent with the correlation func- 
tion determined by Vikhlinin & Forman (1995) for sources 
identifled within ROSAT PSPC deep pointings. A direct 
comparison between the Chandra result with that of Eq. 
(6-3) is complicated by the more than one hundred times 
difference in the angular scales of the two analyses. It is 
doubtful that a single powerlaw model is adequate over 
this range. Furthermore, one is a luminosity ACF while 
the other is a flux limited, number count ACF. Relating 
the two requires understanding the luminosity function 
and its evolution as well as how the X-ray bias depends on 
scale. For these reasons, a direct comparison would be dif- 
ficult to interpret. We only note in passing that the small 
angular scale ACF is a factor of eight larger than that of 
Eq. (6-3) assuming a dependence. 

Finally, the recent harmonic analysis of the HEAOl A2 
data by Scharf et al. (2000) yielded a positive detection 
of structure in the XRB out to harmonic order I ~ 10. 



The present analysis looks at similar maps, so the results 
should be comparable. A direct comparison is complicated 
by the differences in analysis techniques, masking and cor- 
rections to the map. A rough comparison can be made by 
performing a Legendre transform on the ACF model 
of Eq. (6-3). The ACF can be expressed in terms of Leg- 
endre polynomials as 



"(«) 



l)CiPi{cos9) 



(6-4) 



where the Ci constitute the angular power spectrum. Tak- 
ing the Legendre transform 

»i 



Cf = 27r 



uj{9) Pi{cos9) d{cos9) 



(6-5) 



where Pe{9) is the Legendre polynomial of order £. Substi- 
tuting the 9^^ model into this expression results in power 
spectrum coefficients, ~ 4 x 10~^/^ for £ ~ 5. Note 
that for £ 5^ this expression is relatively insensitive to 
the index e in the expression for ujmtr- While these values 
are highly uncertain, they are comparable to those found 
by Scharf et al. (2000) when the sky coverage and differ- 
ences in notation are accounted for. Considering the low 
signal to noise of the data as well as the differences in the 
two analyses, a more detailed comparison would not be 
particularly useful. 

6.2. The Dipole 

Scharf et al. (2000) also searched for the XRB dipole 
using a HEAOl A2 map and similar methods to those 
described in §4. They claim a detection of an intrinsic 
dipole with amplitude A ~ 0.0065, though with a rather 
large region of uncertainty, i.e., 0.0023 < A < 0.0085, 
and in a direction about 80° from that of the Compton- 
Getting dipole, in the general direction of the Galactic cen- 
ter. They used the 3° x 1.5° HEAOl A2 map restricted 
to regions further than 22° from the Galactic plane. In 
addition, regions about sources with fluxes greater than 
2 X 10~^^erg s~^cm~'^ were cut from the map. The sky 
coverage and the level of source removal closely correspond 
to those of our Map C. However, since we used a combi- 
nation of the 3° X 1.5° and 3° x 3° maps, our map has 
significantly less (^ l/\/3) photon shot noise. While our 
analyses are similar, there are some significant differences: 
they corrected the map beforehand for linear instrument 
drift and Galaxy emission while we fit for those compo- 
nents simultaneously with the dipole and with emission 
from the local supercluster which they ignore. 

The upper limit on the intrinsic dipole we find is about 
the same amplitude as the Compton-Getting dipole, i.e., 
A < 0.0052 at the 95 % C.L. Thus we exclude roughly the 
upper half of the Scharf et al. range and believe that their 
claim of a detection is probably an overstatement. While 
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Scharf et al. do not take into account emission from the 
plane of the local supercluster, even if we leave that com- 
ponent out of our fit, the dipole moment (including the 
C-G dipole) increases by only 0.005 TOT cts and is 
still consistent with the C-G dipole. The upper limit on 
the intrinsic dipole with this fit is determined primarily by 
the noise in the fit and is not significantly different from 
that value given above. It is difficult to understand how 
their quoted errors could be two to four times less than 
those quoted in Table |3[ The shot noise variance in the 
map they used was three times greater than in our combi- 
nation map and so one would expect their errors would be 
somewhat larger than those above. They performed only 
a four parameter fit (offset plus dipole) which would result 
in a slight reduction of error; however, this is a bit mis- 
leading since their Galaxy model and linear time drift are 
derived from essentially the same data set. It is possible 
that their lower errors could result from ignoring correla- 
tions in the noise (our detected ACF). In any case, we find 
no evidence for an intrinsic dipole moment in the XRB. 

7. IMPLICATIONS FOR COSMOLOGY 
7.1. Clustering and Bias in the X-ray Background 

The observed X-ray auto-correlation can be compared 
to the matter auto-correlation predicted by a given cos- 
mological model. The linear bias factor for the X-rays can 
then be determined by normalizing to the observed CMB 
anisotropies. Since X-rays arise at such high redshifts, the 
fluctuations we measure are on scales A ^ lOO/i^^Mpc, 
comparable to those constrained by the CMB, i.e., on 
wavelengths that entered the horizon about the time of 
matter domination. 

The predicted X-ray ACF depends on both the cosmo- 
logical model and on the model for how the X-ray sources 
are distributed in redshift, which is constrained by ob- 
served number counts and the redshift measurements of 
discrete sources. We use the redshift distribution de- 
scribed in Boughn, Crittenden and Turok (1998), based 
on the unified AGN model of Comastri et al. (1995). (See 
also the more recent analysis by Gilli et al. (2001).) While 
we will not reproduce those calculations here, the basic 
result is that the XRB intensity is thought to arise fairly 
uniformly in redshift out to z = 4. Our results here are 
not very sensitive to the precise details of this distribution. 

Another issue in the calculation of the power spectrum 
is the possible time dependence of the linear bias. Some 
recent studies indicate that the bias is tied to the growth 
of fluctuations and may have been higher at large redshift 
(Fry 1996, Tegmark & Peebles 1998). For the purposes of 
the power spectrum, an evolving bias will have the same 
effect as changing the source redshift distribution. Again, 
our results are not strongly dependent on these uncertain- 
ties, but they comprise an important challenge to using 
the X-ray fluctuation studies to make precision tests of 
cosmology. 

Figure shows the predicted XRB power spectrum, nor- 
malized to our observations. On the scales of interest, the 
predicted spectra are fairly featureless, and reasonably de- 
scribed by a power law in Ci oc which corresponds 
to the correlation of the form a'(0) = {Oo/Oy. For the mod- 
els of interest, l.I < e < 1.6 for £ < 100, and decreases 
at higher i (smaller separations). Note that the spectra 



calculated by Treyer et al. (1998) appear to be consistent 
with our findings, suggesting e — 1.2. The precise index e 
depends on the position of the power spectrum peak which 
is determined by the shape parameter, F ~ i^mh. Larger 
values of F imply more small scale power and thus higher 
e. 

For simplicity, we normalize to the X-ray correlation 
function at 4.5°, w(4.5°) = 1.0 ± 0.25 x lO"'^. This sepa- 
ration is large enough to be independent of the PSF con- 
tribution to the ACF, but not so large that the attenu- 
ation from the large scale fits becomes significant. Also, 
the value of the fit ACF at 4.5° is nearly independent of 
the index e (see §5). As can be seen from Figure |^, this 
normalization fixes the power spectrum at ^ ~ 5 — 7. 

We normalize the fluctuations to the COBE power spec- 
trum as determined by Bond, Jaffe & Knox (1998). How- 
ever, it should be noted that fits to smaller angular CMB 
fluctuations indicate that using COBE alone may some- 
what overestimate the matter fluctuation level (Lahav et 
al. 2002). The biases derived from the models appear to 
be largely insensitive to the matter density. This is due to 
a cancellation of two effects: the CMB normalization and 
the power spectrum shape (White & Bunn 1995). The bi- 
ases are roughly inversely proportional to h. Typical biases 
appear to be bx — 2.3 ± 0.3(0.7//i)*'-^, increasing shghtly 
as F decreases and the peak of the power spectrum moves 
to larger scales. 

7.2. The Intrinsic X-ray Dipole 

The theoretical models normalized to our observations 
predict the intrinsic power on a wide range of scales, as- 
suming the X-ray bias is scale independent. In particular, 
these models give a prediction for the variance of the in- 
trinsic dipole moment. We can compare our model predic- 
tions to the the upper limit for the intrinsic dipole to see 
if we should have observed it in the X-ray map. 

The dipole amplitude in the z direction is related to the 
spherical harmonic amplitude by — y^S/inaio. Thus, 
the expected dipole amplitude is related to the power spec- 
trum by 

(A2)^3x = Aci. (7-1) 

Note that there is considerable cosmic variance on this, 
as it is estimated with only three independent numbers; 
SCi/Ci = (2/3)"'^/^, which corresponds to a 40 % uncer- 
tainty in the amplitude of the dipole. 

Also shown in Figure ^ is the level of our dipole limit, 
translated using equation (7-1), which corresponds to 
Ci < 3.8 X 10^^. While our limit on the dipole hmit 
is at the 95 % C.L., this translates to a 80 % C. L. limit 
on the variance Ci when cosmic variance is included. As 
discussed above, the 95 % upper limit is four times weaker 
when cosmic variance is included, Ci < 1.5 x 10"''. Nor- 
malized to the ACF, all the theories are easily compatible 
with the Ci bound. The large cosmic variance associated 
with the dipole makes it difficult to rule out any cosmo- 
logical models. 

With our detected level of clustering, typical theories 
would predict a dipole amplitude of A ~ 0.003. While the 
theories are not in conflict with the dipole range claimed 
by Scharf et al., they strongly prefer the lower end of their 
range, even for the most shallow of the models. A dipole 
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amplitude A > 0.005 would be very unlikely from the 
models, indicating either a significantly higher bias than 
we find or a model with more large scale power (e < 1.1). 

7.3. The dipole and bulk motions 

The dipole of the X-ray background provides another in- 
dependent test of the large scale X-ray bias through its re- 
lation to our peculiar velocity (e.g. see Scharf et al. (2000) 
and references therein.) Like the gravitational force, the 
flux of a nearby source drops off as an inverse-square law, 
so the dipole in the X-ray flux is proportional to the X- 
ray bias times the gravitational force produced by nearby 
matter. Our peculiar motion is a result of this force, and is 
related to the gravitational acceleration by a factor which 
depends on the matter density. 

In typical CDM cosmologies, the dipole and our pecu- 
liar velocity arise due to matter at fairly low redshifts 
(z < 0.1). If this is the case, it is straight forward to 
relate their amplitudes. Following the notation of Scharf 
et al, we define = / df^/(n)n" = 47r/A°/3. Using lin- 
ear perturbation theory, one can show the local bulk flow 
is 

where px (0) is the local X-ray luminosity density, bx (0) is 
the local X-ray bias and / ~ fl^f is related to the growth 
of linear perturbations (Peebles 1993). From the mean 
observed intensity (Gendreau et al. 1995) and the local 
X-ray luminosity density (Miyaji et al. 1994) we find that 
/ ~ 2.4px(0)c/47ri7o. This implies that 

|?;| ~2.4x lO^A-^^kms-i. (7-3) 
ox[0) 

This relation was derived by Scharf et al. (2000), though 
their numerical factor was computed from a fiducial model 
rather than directly from the observations, as above. 
In any case, the uncertainty in px(0) is considerable, 
6 X lO^'^ er g s~^Mpc'^ < px{0) < 15 x 10^^ erg s"'^ Mpc~^ 
(Miyaji et al. 1994), and so is the uncertainty in this re- 
lation. 

Our maps have bright sources removed, which corre- 
spond to nearby sources out to 60h~^Mpc. Thus, we need 
to compare our dipole limit to the motion of a sphere of 
this radius centered on us. Typical velocity measurements 
on this scale find a bulk velocity of ugo — 300 ± 100 km/s 
(see Scharf et al. (2000) for a summary.) With our dipole 
limit, this implies that fl'^^/bx{0) > 0.24±0.08 where the 
uncertainty in l/px{0) is not included. This constraint is 
independent of cosmic variance issues. While the diam- 
eter of the local (Virgo) supercluster is generally consid- 
ered to be in on the order of AOh~^ to 50h~^ Mpc (e.g., 
Davis et al. 1980), there is evidence that the overdensity 
in the Supergalactic plane extends significantly beyond 
60/i~i Mpc (Lahav et al. 2000). One might, therefore, 
suspect that our correction for emission from the local su- 
percluster might effectively remove sources at distances 
greater than 60h~^ Mpc in that plane. In any case, that 
correction made very little difference in dipole fits (see 
§6.2) so our conclusions remain the same. 

Note that this limit could potentially conflict with 
previous determinations by Miyaji (1994) who found 
^^/bx{0) — /45/3.5, where /45 is the fraction of grav- 
itational acceleration arising from R < A5h~^ Mpc. This 



is consistent only for f^^ ^ 1, which is larger than is usu- 
ally assumed (/45 ^ 0.5). However, this limit comes from 
studies of a fairly small sample (16) of X-ray selected AGN 
and is subject to signiflcant uncertainties of its own. 

For typical biases suggested by the observed clustering 
{bx ~ 2.3), our constraint suggests a somewhat high mat- 
ter density, fi™ > 0.37, for wgo — 300 km/s. This is consis- 
tent with the ISW constraint discussed below and also with 
previous analyses of bulk velocities which tend to indicate 
higher flm- However, if the bulk velocity is smaller and/or 
px(0) larger, this constraint is weakened. In addition, we 
have assumed a constant X-ray bias. If the bias evolves 
with redshift, then the local value could be considerably 
smaller which would also weaken this bound. 

7.4. The Integrated Sachs- Wolfe Effect and Q\ 

In models where the matter density is less than unity, 
microwave background fluctuations can be created very 
recently by the evolution of the linear gravitational poten- 
tial. This is known as the late time integrated Sachs- Wolfe 
(ISW) effect. Photons gain energy as they fall into a po- 
tential well, and loose a similar amount of energy as they 
exit. However, if the potential evolves significantly as the 
photon passes through, the energy of the photons will be 
changed, leaving an imprint on the CMB sky. The spec- 
trum is modified most on large scales where the photons 
receive the largest changes. 

The CMB anisotropics created in this way are naturally 
correlated with the gravitational potential. Thus, we ex- 
pect to see correlations between the CMB and tracers of 
the local (z ~ 2) gravitational potential such as the X-ray 
background (Crittenden & Turok 1996). These correla- 
tions are primarily on large scales such as those probed by 
the HEAO survey. 

In an earlier paper, we searched for a correlation be- 
tween the HEAO maps and maps of the CMB sky pro- 
duced by COBE. We failed to find such a cross correlation 
and were able to use our limit to constrain the matter 
density and the X-ray bias (Boughn, Crittenden & Turok 
1998, hereafter BCT). However, translating our measure- 
ment into a cosmological bound was ambiguous because 
the level of the intrinsic structure of the XRB was un- 
known at the time. With the observation of the X-ray 
ACF presented here, we are in a position to revisit the 
cosmological limits implied by these measurements. 

To make cosmological constraints, we compare the ob- 
served X-ray /CMB cross correlation to those predicted by 
KCDM models. As above, we normalize the CMB fiuc- 
tuations using the band powers of COBE (Bond, Jaffe & 
Knox 1998) and also normalize the X-ray fluctuations as 
discussed in §7.1. The cross correlation analysis of BCT 
was performed with a coarser pixelization (2.6° x 2.6°) 
than the ACF discussed above. We include this effect by 
using the numerically calculated pixelization window func- 
tion. The COBE PSF used was that found by Kneissl & 
Smoot (1993) and we used a 2.9° FWHM Gaussian for the 
underlying X-ray PSF (recaU that the 3.04° FWHM beam 
found above includes a 1.3° x 1.3° pixelization.) 

The calculation of the HEAO-COBE cross correlation 
was discussed in BCT and has not changed. The results 
are shown in Figure ^, along with predictions for three dif- 
ferent values of f^A- While the X-ray bias depends strongly 
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on the Hubble parameter, the predicted cross correlation 
is only weakly dependent on it, changing only 10% for 
reasonable values of Hq. The cross correlation depends 
primarily on no correlation is expected if there is no 
cosmological constant and the ISW effect increases as Ha 
grows. The error bars in Figure ^ are calculated from 
Monte Carlo simulations and arise primarily due to cos- 
mic variance in the observed correlation. The error bars 
are significantly correlated. 

The observed correlation is most consistent with there 
being no intrinsic cross correlation (f^A = 0.0). We set 
limits by calculating the likelihood of a model relative to 
this no correlation model. Using the frequentist criterion 
used in BCT, flj^ < 0.65 at the 98% C.L., f^A < 0.60 at the 
95%. C.L. Almost identical limits arise from a Bayesian 
approach, where the relative likelihoods are marginalized 
over, assuming a constant prior for fl^ > 0. Figure || 
shows a one-dimensional slice through the likelihood sur- 
face, where only the cross correlation information has been 
used to calculate the likelihood. 

One of the major assumptions we made in interpreting 
the above result is how the sources of the XRB are dis- 
tributed in redshift. It is likely that current models of the 
luminosity function will have to be substantially modified 
as further deep observations of the sources of the XRB 
are made. However, as pointed out above, the ISW is 
relatively insensitive to the exact shape of the redshift dis- 
tribution of luminosity. If the true distribution includes a 
substantial fraction of the luminosity at redshifts greater 
than 1, then the above results will not change dramatically. 
On the other hand, our constraint on Q\ is quite sensitive 
to the value of the bias parameter. If the sources of the 
XRB should turn out to be unbiased, i.e., 6x = 1, then 
the constraint on f^A could be weakened dramatically. We 
hasten to add that such a low bias would require that the 
ACF of Figure 5 be reduced by more than a factor of four, 
which seems unlikely. Previous determinations of X-ray 
bias have resulted in a wide range of values, 1 < 6a: < 7 
(see Barcons et al. 2000 and references therein). Its clear 
that firming up the value of bx and determining how it 
varies with scale and redshift will be required before the 
ISW effect can be unambiguously interpreted. 

The above limit may be compared to what we found 
from cross correlating COBE with the NVSS radio galaxy 
survey (Boughn & Crittenden 2002). There we also found 
no evidence for correlations, and were able to put a 95 
% C.L. Hmit of l^A < 0.74, with some weak dependence 
on the Hubble constant. While the above limit provides 
important confirmation of that result, it should be noted 
that these two limits are not entirely independent. Radio 
galaxies and the X-ray background are, indeed, correlated 
with each other (Boughn 1998). 

An important source of noise in the cross correlation of 
Figure 7 is instrument noise of the COBE DMR receivers. 
In addition, the relatively poor angular resolution of the 
COBE radiometers reduce, somewhat, the amplitude of 
the ISW signal. Therefore, some improvement can be ex- 
pected by repeating the analysis on future CMB maps, 
such as that soon to be produced by NASA's MAP satel- 
lite mission. If such an analysis still finds the absence of 
an ISW effect, then the current AC DM model would be 



in serious conflict with observational data if the X-ray bias 
can be similarly constrained. On the other hand, a posi- 
tive detection would provide important evidence about the 
dynamics of the universe even if the X-ray bias remains 
uncertain. 

8. CONCLUSIONS 

By carefully reconstructing the HEAO beam and 
analysing its auto-correlation function, we have been able 
to confirm the presence of intrinsic clustering in the X-ray 
background. This gives independent verification of the 
multipolc analysis of Scharf et al. (2000) and the level of 
clustering we see is comparable. The clustering we see is 
in excess of that predicted by standard cold dark matter 
models and indicates that some biasing is needed. The 
amount of biasing required depends on the cosmological 
model and on how the bias evolves over time; if the bias is 
constant, typical models indicate that 5x — 2. The biases 
of galaxies, clusters of galaxies, radio sources, and quasars 
have yet to be adequately characterized and so whether or 
not the above X-ray bias is excessive is a question that, 
for the present, remains unanswered. 

We have also confirmed, at the 2-3 a level, the detec- 
tion of the Compton-Getting dipole in the X-ray back- 
ground due to the Earth's motion with respect to the rest 
frame of the CMB. However, we have been unable to con- 
firm the presence of an intrinsic dipole in the XRB and 
have actually been able to exclude a significant part of 
the range reported by Scharf et al. (2000). While our 
dipole limit is still too small to confiict with any of the 
favored CDM models, combining our dipole limit with ob- 
servations of the local bulk flow enable us to constrain 
fl^/bx{0) > 0.24. For constant bias models, this sug- 
gests a relatively large matter density, as is also seen in 
for other velocity studies; however, the uncertainty in this 
limit is still considerable. 

With the observed X-ray clustering, large A — CDM 
models predict a detectable correlation with the cosmic 
microwave background arising via the integrated Sachs- 
Wolfe effect. That we have not observed this effect sug- 
gests f^A ^ 0.60. This is beginning to conflict with models 
preferred by a combination of CMB, LSS and SNIA data 
(e.g., de Bernardis et al. 2000 & BahcaU et al. 1999). 

This work gives strong motivation for further obser- 
vations of the large scale structure of the hard X-ray 
background. Better measurements of the full sky XRB 
anisotropy are needed, as is more information about the 
redshift distribution of the X-ray sources. This will be es- 
sential for cross correlation with the new CMB data from 
the MAP satellite and to bridge the gap between the CMB 
scales and those probed by galaxy surveys such as 2-dF 
and SDSS. 
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Fig. 1. — The combined map from the HEAOl A2 medium and high energy detectors, pixelized using the standard COBE quad cubed 
format (1.3° X 1.3° pixels.) The effective beam size is approximately 3°. The most visible features, the Galactic plane and the nearby bright 
sources, are removed from the maps we analyze. 
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Fig. 2. — The mean point spread function for the combined map found by averaging the individual PSFs of sixty strong HEAOl point 
sources. The data is well fit by a Gaussian with FWHM of 3.04°. 
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Fig. 3. — The auto-correlation function of the HEAOl A2 map with bright sources and the Galactic plane removed and corrected for 
large-scale, high Galactic latitude structure. The dashed curve is that expected from beam smearing due to the PSF of the map while the 
solid curve includes a contribution due to clustering in the XRB (see §5). 
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Fig. 4. — The residuals of the ACF fit from Figure 3, after the shot noise, PSF and a simple model of the intrinsic fluctuations have been 
removed. 




6 ( degrees ) 

Fig. 5. — The intrinsic ACF, with shot noise and PSF fits removed. For comparison, a simple model for the intrinsic correlations is 
shown. The data beyond 9° is not used because of uncertainty due to the fitting of the large scale structures. The model has been smoothed 
by the PSF and corrected for the removal of the large scale structures, which suppresses the correlations on scales larger than 10°. 
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Fig. 6. — The power spectrum for a range of cosmologies normalized to the observations (Hq = 70 km Mpc~^ .) The various cosmologies 
show a range of slopes, from 1.1 < e < 1.6 and the observations fix them at £ ~ 5. Also shown is the 95 % upper limit from the dipole, 
excluding cosmic variance. With cosmic variance, the limit shown is at the 80 % confidence level, and the 95 % upper limit is four times 
higher. The green line shows the suppression arising from beam smoothing, which smoothes scales £ > 50. 
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Fig. 7. — The calculated X-ray/CMB cross correlation. The error bars are highly correlated. Also shown are the predictions for three 
A - CDM models with varying Qa (^^^0 = 70 km s-^Mpc-^.) 
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Fig. 8. — The relative probability of the observed cross correlation for varying cosmological constant, with the Hubble constant fixed 
(Ho = 70 km Mpc~^ .) The best fit is for no correlation. 



